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ABSTRACT 

We present a detailed, self-consistent model of radia- 
tively driven stellar outflows which couples the radia- 
tive transfer and hydrodynamics equations. The circum- 
stellar envelope, which consists of gas and dust, is de- 
scribed as a two-component fluid to account for rela- 
tive drifts. Our results agree with both molecular line 
observations and infrared continuum spectra, and show 
that steady-state outflows driven by radiation pressure 
on dust grains adequately describe the surroundings of 
the majority of cool luminous evolved stars. 

1. INTRODUCTION 

Astronomical objects are often characterized by 
large linear scales, huge masses and powerful flows of 
energy. Such environments emphasize different physical 
effects than those encountered in a laboratory or in en- 
gineering applications. An example is the radiation pres- 
sure force. This force results from the transfer of momen- 
tum from the radiation field to the intervening medium. 
For isotropic interaction between the radiation and the 
matter, the acceleration due to radiation pressure is ([1]) 
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where c is the speed of light, Fa is radiative flux vector, 
X\ is the opacity per unit mass of the medium (includ- 
ing both absorption and scattering) and A is wavelength. 
In the engineering applications this force is usually neg- 
ligible, but in some astronomical environments it can 
dominate over all other forces. 

Cool luminous evolved stars are one astronomical 
environment where radiation pressure plays a major 
role. These stars, which are in the last stages of their 
evolution, experience mass- loss. Although the exact na- 
ture of the mechanism that governs the initial ejection 



of mass from the stellar surface is still not clear [2] , sub- 
sequent outflow is driven by the radiation pressure force 
[3,4]. When the expanding gas reaches a certain distance 
from the star it becomes sufficiently cool that dust grains 
(solid particles with sizes ^ 1/im) begin to condense. 
Radiation pressure on the dust continues to push the 
material away from the star and a circumstellar shell, or 
envelope, is formed. Since the dust grains absorb stel- 
lar radiation and reradiate it in the infrared, the overall 
radiative spectral distribution is shifted toward longer 
wavelengths and many of these objects are invisible in 
the optical part of the spectrum. 

In order to solve the radiation transfer in such 
envelopes, the mass density distribution needs to be 
known. However, a given density distribution implies a 
particular velocity field which, in turn, is determined by 
the radiation pressure force. Thus, the radiative transfer 
and the dynamics are inherently coupled and have to be 
solved simultaneously. 

Our model assumes spherical symmetry and cou- 
ples the momentum conservation and radiative transfer 
equations. The circumstellar envelope, which consists of 
gas and dust, is treated as a two-component fluid. We 
present an investigation of the dynamics and spectra of 
dusty envelopes around late-type stars, and discuss the 
implications for observations. 

The rest of this paper is organized as follows: first 
we discuss the conditions for the radiation pressure force 
to be important (§2), then we describe our model (§3), 
and present the results (§4). Some limitations of the 
model are discussed in §5. 

2. WHEN IS THE RADIATION 
PRESSURE FORCE IMPORTANT? 

Eq. (1) shows that the acceleration due to the radi- 
ation pressure, a rp , is proportional to the radiation flux. 
This flux cannot be arbitrarily large because dust grains 
are radiatively heated, and in sufficiently intense radia- 
tion field can evaporate. Thus the upper bound on the 
radiation pressure force is set not by the intensity of the 
radiation field but by the evaporation temperature of 
the grains, T cvap -1000 K. 

The radiation pressure force is maximized when the 
grains are heated by a directional radiation field, |Fa| = 
c u\, where u\ is the radiation energy density. In this 
case, the total (bolometric) flux absorbed by grains when 
they are heated to T cvap is: 

In deriving the numerical constant we have assumed 
small grains with opacity proportional to A -1 . With this 
flux, the acceleration due to the radiation pressure is 

max ( T cva p \ _2 ,„* 
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where g = 9.81 m s -2 . This is the maximal acceleration, 



that a grain can experience due to the radiation pressure 
force, and is independent of the radiation field intensity. 

Dust grains are usually mixed with gas particles. In 
such a mixture the momentum which grains gain from 
the radiation field is transferred to gas particles via col- 
lisions [5]. In this case the effective acceleration of such 
a two-fluid system has to be multiplied by the dust-to- 
gas mass ratio (<C 1), and thus a rp <C g. As a result, in 
engineering applications the radiation pressure force is 
always negligible compared to the gravitation force. 

When is the radiation pressure force important in 
astrophysical systems? The radiation pressure force will 
drive an outflow when it overcomes the gravitation force. 
Taking a typical circumstellar value for the dust-to-gas 
mass ratio of 0.002 [6], a rp ~10~ 4 g. Furthermore, this 
value can be even 10-100 times smaller when the dust 
drift effects become significant (this correction is propor- 
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tional to p gaB , and puts a lower limit on the mass-loss 
rate [2], see §3.1). 

The gravitational acceleration is 
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where G is the gravitational constant, M is the stellar 
mass, and r is the distance from the star. Dust grains 
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where L* is the stellar luminosity, and Lq= 3.86 x 10 
W m -2 is the solar luminosity. The gravitational accel- 
eration at the dust evaporation radius is thus 
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where M = M*/M Q (M©= 2 x 10 30 kg is the solar 
mass), and L — L^/Lq. Note that both a rp and a g scale 

With Tetap. 

As the stellar luminosity increases the dust evapo- 
ration radius is pushed further out, and a g decreases. 
For sufficiently small M/L, a g becomes smaller than 
a rp , and the radiation pressure force drives an outflow of 
dust and gas. Taking a typical value of M~l, the mini- 
mum required luminosity is 10 3 to 10 4 L@. This is why 
a significant luminosity increase during the late stages 
of stellar evolution is accompanied by outflows; it is the 
small local a g , rather than the large a rp , which enables 
the radiatively driven stellar outflows. 



3. A MODEL FOR RADIATION 
DRIVEN OUTFLOW 

3.1 The Mathematical Model 

In this section we briefly present the basic equations 
describing the mathematical model. Detailed discussion 
of the model and the simplifications can be found in 
[2,7]- 

The equation of motion for the gas is 
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where v is the gas velocity. The first term on the right- 
hand side of eq. ([?]) describes the radiation pressure 
force, the second term is the gravitational force exerted 
by the star. The force due to the gas pressure gradient 
is omitted since it is negligible for the highly supersonic 
outflows described here. In cool stars, the radiation pres- 
sure force acts only on the dust grains (gas opacity is 
negligible). Since the dust flows out faster than the gas, 
the dust opacity decreases in proportion to the ratio of 
the gas and dust outflow velocities 
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where the dust drift velocity with respect to the gas flow 
is determined from 
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Here Qf is the flux averaged absorption efficiency, and 
M is the mass-loss rate (for a detailed discussion of dust 
drift see [7]). 

At first it appears from eq. (Q) that outflows with 
arbitrary mass-loss rates can be driven by the radiation 
pressure force. However, the allowed mass-loss rate is 
limited from both below and above. The lower limit is 
set by the above mentioned drift effects: as the mass-loss 
rate, i.e. gas density, decreases, the dust-gas coupling de- 
creases and with it the effective opacity. In the limit of 
large mass-loss rates, practically the entire stellar radia- 
tion is absorbed by dust, and shifted to longer, infrared 
wavelengths. As opacity is in general a decreasing func- 
tion of wavelength, the effective opacity decreases with 
mass-loss rate in this regime, resulting in an upper limit 
on the mass-loss rate. 

3.2 The Solution Method 

The coupled system of the dynamics and radiative 
transfer equations is solved by an iterative scheme. For a 
given velocity field, the radiative transfer is solved sepa- 
rately by an exact scheme described in [8]. With known 
radiation flux, an updated velocity field is calculated 
from eq. (Q), and these steps are repeated until both the 
radiation and velocity fields converge. As an initial ve- 
locity field we take an analytic approximation obtained 
by neglecting drift and reddening effects, and gravita- 
tional force [9]. 

The radiative characteristics of the star are largely 
irrelevant for the solution. We assume that the star ra- 
diates as a black body with T« = 2500 K, typical for 
late-type red giant stars. Similarly, we assume that the 
dust evaporation temperature T eV ap= 1000 K. The dust 
optical properties are important because they determine 
the detailed signature of the emitted spectrum. In this 
work we consider the two major grain types that have 
been proposed for late-type stars. Dielectric coefficients 
for "astronomical silicate" are taken from Ossenkopf et 
al. [10], and for amorphous carbon from Hanner [11]. 
A more detailed study including 5 different grain types 
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Figure 1. Model results for the dependence of outflow veloc- 
ity on mass-loss rate M, luminosity L*=Z/4l0 4 Lq, and gas- 
to-dust ratio r gt j. The plotted variables, whose forms follow 
from the model scaling properties and dimensional analysis, 
describe families of self-similar solutions for two types of dust 
grains, as marked. 



and various mixtures is presented in [4]. Absorption 
and scattering efficiencies are evaluated using Mie the- 
ory, assuming the dust grains are homogeneous spheres 
following the size distribution proposed by Mathis et 
al. [12]. Carbon grains have smooth opacity, decreasing 
with wavelength, while astronomical silicate has charac- 
teristic peaks at 9.7 /im and 18 /im. 



4. RESULTS 

4.1 Scaling Properties of the Solution 

For given dust optical properties, stellar tempera- 
ture, dust evaporation temperature, gas-to-dust ratio, 
mass-loss rate, and luminosity, our model determines 
the outflow velocity field. Optical properties and the two 
temperatures are usually well constrained, so in the re- 
mainder of this paper we will consider gas-to-dust ratio 
r gd , v e , M, and L* as the relevant quantities describing 
a model. 

The velocity field can be described by terminal out- 
flow velocity v e , which sets the scale, and shape u — 
v/v e . From the scaling properties of the system, v e can 
be expressed as the product of a velocity scale propor- 
tional to L t 1/ ^ 4 r gd 1 ^ 2 , and a dimensionless function of the 
dust optical depth [9] . Because the solution is character- 
ized by a single parameter, it follows from dimensional 
analysis that the terminal outflow velocity is related to 
the input quantities through some (quite involved) func- 
tion $ such that 



1/4-1/2 



$(- 



M 



* ' gd 



) 



(10) 



The argument of function 3> follows from the dependence 
of optical depth on M, L t , and r gd . That is, for given 




Figure 2. Model results for the shape of velocity field u = 
v(r)/v e for two types of dust grains, as marked. Different 
curves correspond to different dust optical depths at 0.55 fim 
(dotted lincs:l, dot-dashed lines: 10, dashed lines: 100). The 
solid line corresponds to an approximate solution given by 
cq. (11) with k = 2. 



grains, there is a single family of self-similar solutions, $, 
which relates v e r^L, ~ 1/4 to ML, ~ 3/4 r , 1/2 . Strictly 



speaking, 3> also depends on M* , and e = vi /v e , where 
vi is the velocity at the inner boundary, ri , but detailed 
model calculations show that this dependence is negli- 
gible for all practical purposes. Solutions for two types 
of grains considered in this work are shown in figure 1 
(obtained for e = 0.1 and M» = 1 M ). These solutions 
can be used to reconstruct the v e (M) relationship for 
arbitrary L» and r gd . 

We find from detailed models that the shape u is 
well described by 



«=(e fc + (l-e fc )(l-^)) 



l/k 



(11) 



The characteristic power index k is a weak function of 
dust optical depth (fig. 2), ranging from 1.5 in optically 
thin envelopes to ~2.5 in optically thick ones (k=2 cor- 
responds to an idealized case with acceleration propor- 
tional to r~ 2 ). Again, the dependence on M* and e is 
negligible. 
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Figure 3. Comparison of v e and M obtained from molecular 
line observations (symbols) with model predictions (lines). 
For more details please see 4.2. 



4.2 Dynamical Properties of the Solution 

The predicted relationship between the terminal 
outflow velocity v e , and mass- loss rate M, can be used 
to compare the model to observations. Observations of 
molecular line emission (both thermal and maser) di- 
rectly measure v e , and with the aid of suitable emis- 
sion models constrain M. Straightforward comparison of 
model predictions with observed v e and M is hampered 
by uncertain values of L* and r g d- Consequently, a direct 
utilization of the results presented in figure 1 is imprac- 
tical because effects of changing L t and r g d are coupled 
together. Their impact can be disentangled by compar- 
ing the observed and predicted values of Mi) e " r~ d and 
MveL,^ 1 . Such a diagram is shown in figure 3. Lines are 
model predictions and symbols are the observed values 
for two samples of stars: oxygen rich (O) stars which are 
associated with silicate dust, and carbon rich (C) stars 
associated with carbon dust. We have assumed that all 
stars have L4 = L,/lO 4 I/0 = 1, and constrained r g d 
by requiring best agreement (in a least-squares sense for 
the whole sample) between the modeled and observed 
values of Mv^~ 3 r~ d ■ Best values are 250 for O stars, and 
400 for C stars. It is evident from figure 3 that these 
values are practically independent of assumed value for 
luminosity because different values of L* shift the data 
points only horizontally. 

The obtained values for r g d are in agreement with 
other independent estimates (e.g. [6]). Close agreement 
between the model predictions and observations pro- 




Figure 4. Model spectra (normalized \F\) for astronomical 
silicate grains for a sequence of flux averaged optical depths 
Tp. Solid lines present the spectral shape of the emergent 
radiation, the sum of two components — the dust emission, 
indicated in dotted lines, and the emerging stellar radiation, 
the dashed lines. The model prediction for the IRAS [25]- 
[12] color is indicated. In each panel, the inset displays the 
measured spectrum of an actual source whose IRAS number 
is marked above the top left corner. The measured [25]— [12] 
color of the source is listed in the inset in large type. 



vides strong evidence that outflows around late-type 
stars are radiatively driven. Additional evidence is pro- 
vided by the analysis of dust infrared emission. 

4.3 Infrared Emission 

An important property of the spherically symmet- 
ric radiative transfer problem is that, for given grains, 
the dimensionless spectral shape fx — (\F\)/F, where 
F = J °° F\d\ is the bolometric flux, depends only on 
the radial density profile and overall optical depth [8]. 
The radial density profile for an outflow is fully deter- 
mined by the velocity profile through the mass conserva- 
tion relation. Since the velocity profile depends only on 
optical depth, it follows that for each grain type there is 
a family of spectra parametrized by optical depth. That 
is, irrespective of their individual values, any combina- 
tion of relevant parameters produces the same spectral 
shape as long as the corresponding optical depth does 
not change. 

A detailed comparison of observations with model 
spectra is given in [4], and here we only present results 
for silicate grains as an illustration. The six panels of 



figure 4 display a series of spectra parametrized by the 
flux averaged optical depth 

T F = J T X ^d\. (12) 

The plotted wavelength range matches the spec- 
tral observations obtained by the Infrared Astronomical 
Satellite Low Resolution Spectrometer (IRAS LRS) [13]. 
Dashed lines represent the direct stellar radiation emerg- 
ing from the system, dotted lines are the dust emission 
and full lines are the total emission — the model pre- 
dictions of observed spectra. In optically thin envelopes 
{tf < 1) the strength of the 9.7 /jm silicate emission 
feature increases with optical depth until tf reaches ~ 
1. With further increase in tf, this emission turns into 
an absorption feature whose depth is proportional to 
tf- As an illustration, the inset in each panel displays 
the observed spectrum of an actual source whose IRAS 
number is indicated above the top left corner. Each panel 
also indicates the [25]-[12] coloijj of the computed model 
while the number listed in large type in the inset is the 
IRAS color of the displayed source. The close agreement 
of spectral evolution between model predictions and ob- 
servations is evident. 

5. DISCUSSION 

Since our models assume steady state, they can only 
describe the time-averaged behavior of the outflows. 
This is an adequate description as long as time variabil- 
ity occurs only on time scales shorter than the averaging 
time. Calculations show that most of the radiative ac- 
celeration takes place within r ~ 10 16 cm. With typical 
outflow velocities of ~ 10 km s _1 , this distance is cov- 
ered in ~ 300 years, the relevant time scale for averaging 
density. Although all late-type stars display pulsational 
or irregular variability, those variations are character- 
ized by time-scales of a few years at most, much shorter 
than the averaging time. Thus the assumption of steady 
state seems justified. Additional time dependence can 
be introduced by chemical evolution of the stellar atmo- 
sphere, possibly leading to gradients in the grain chemi- 
cal composition. Though our models do not include such 
an effect, they adequately describe the envelopes as long 
as chemical gradients occur only at r <; 10 16 cm. The 
reason is that the grain composition in that region is 
largely irrelevant since all grains behave similarly at low 
temperatures and long wavelengths. And since spatial 
gradients with length scale ^ 10 16 cm imply temporal 
variations with time scales 100 years, the durations of 
phases not covered by our models do not exceed ~ 100 



yearsPl. If typical lifetimes in the mass-losing phase are 
~ 10 -10 s years, only ~ 0.1-1% of all sources cannot 
be described by our models. Evolutionary effects char- 
acterized by time scales longer than ~ 100 years can be 
incorporated into our models as adiabatic changes of the 
relevant sources. 

An important ingredient of our models is the as- 
sumption of spherical symmetry. Small departures from 
sphericity, such as a slight elongation, should not affect 
the results significantly. However, our models do not ap- 
ply in the case of major deviations from spherical sym- 
metry such as disk geometry or strong dumpiness with 
characteristic size scale less than ~ 10 16 cm. The suc- 
cess of our models in describing the infrared spectra and 
mean mass-loss rates indicates that such deviations may 
not be important for the majority of late-type stars. 
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* Colors are convenient parameters used in astronomy 
to describe spectral shapes. Here we define [25]-[12] as 
log(i ? 2s/i ; l2), where F%2 and F25 are fluxes at 12 /im and 
25 fim, respectively. 



t It has been recently suggested that mass-loss rate for some 
stars might change on similar time scales [14] . 



